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Abstract. Photon emission by an electron embedded in a strong external field of general form is studied theoretically. 
The external field considered is a plane wave electromagnetic field of any number of components, period and polar- 
isation. Exact, Volkov solutions of the Dirac equation with the 4-potential of the general external field are obtained. 
The photon emission is considered in the usual perturbation theory using the Volkov solutions to represent the elec- 
tron. An expression for the transition probability of this process is obtained after the usual spin and polarisation sums, 
trace calculation and phase space integration. The final transition probability in the general case contains a single sum 
over contributions from external field photons, an integration over one of the phase space components and the Fourier 
transforms of the Volkov phases. The validity of the general expression is established by considering specific external 
fields. Known specific analytic forms of the transition probability are obtained after substitution of the 4-potential for a 
circularly polarised and constant crossed external field. As an example usage of the general result for the transition prob- 
ability, the case of two circularly polarised external fields separated by a phase difference is studied both analytically 
and numerically. 



PACS. 12.20.Ds QED Specific calculations 

1 Introduction 

The High Intensity Compton Scattering (HICS) is the radiation 
of a single photon by an electron in an external electromag- 
netic field. If the external field is sufficiently intense, multi- 
photon interactions between the electron and the external field 
become likely and many terms in the usual perturbation series 
describing the overall interaction, contribute to the final transi- 
tion probability. Consequently, an interaction picture in which 
the interaction with the external field, considered to be classi- 
cal, is calculated exactly and in which the interaction with the 
quantised boson field is perturbative - the Furry picture - is usu- 
ally employed Q]. 

The HICS process is currently the subject of much study 
largely because new experimental facilities are being planned 
which will be able to test specific non-linear features of the the- 
ory [2,3 1 . A variety of experimental conditions can be created 
each of which require a specific analytic form of the HICS tran- 
sition probability in order to study. The aim of this paper is to 
provide a shortcut to this theoretical work by generalising the 
theoretical description of the transition probability of the HICS 
process and to illustrate its specific usage in cases of practical 
interest. 

Much theoretical work on the HICS process, using the Furry 
picture and Volkov solutions, was spurred by the construction 
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of the first LASER in 1960 and the possibility of laboratory 
tests of the theory. The transition amplitude of the HICS pro- 
cess was found to be dependent on the state of polarisation of 
the external field. For the case of a linearly polarised exter- 
nal electromagnetic field, the HICS transition probability con- 
tained an infinite summation of complicated functions which 
were evaluated in limiting cases only |4|. In contrast, a circu- 
larly polarised external field results in a transition probability 
containing Bessel functions, the properties of which are well 
known |5]|6)- A circularly polarised external field introduces 
an azimuthal symmetry into the HICS process which results 
in an analytically less complicated HICS transition probability 
Q. The HICS process for the case of an elliptically polarised 
external electromagnetic field and for the case of two orthogo- 
nal, linearly polarised fields was considered by ISJ. 



An important effect emerging from all this theoretical work 
was a dependency of the energy of the radiated photon (a fre- 
quency shift) on the intensity of the external field EH). The 
existence of this frequency shift mechanism allows the HICS 
process to be used as a generator of high energy photons ifTUl 
[TT1IT21 . The polarisation properties of high energy photons pro- 
duced by the HICS process are of fundamental importance in 
nuclear physics applications, for instance in the study of abnor- 
mal parity components in the deuteron wavefunction [1311141 . 
The polarisation state of the emitted HICS photon, as a function 
of initial polarisation states, for a circularly polarised external 
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field, was studied by B13II15L 

In experimental work on the HICS process, a 1.053 fim 
laser focused to a 5 fj,m spot size to generate a peak laser inten- 
sity of ~ 10 18 Wcm~ 2 interacted with electrons resulting in a 
longitudinal shift in scattered electron momenta due to multi- 
photon contributions from the laser field. Also observed was 
the predicted electron mass shift due to the presence of the ex- 
ternal field B16II17II . The HICS process and another strong field 
process, the one photon pair production process, were studied 
at SLAC in the 1990s by interacting a Nd:glass laser of peak 
intensity 0.5 x 10 18 W cirT 2 with a 46 GeV electron beam [ 18 1 

The prospect in the near future of extreme light sources pro- 
ducing laser pulses with intensities exceeding those of previ- 
ous experiments has prompted new theoretical work in pulsed 
background fields. The HICS process was studied in a circu- 
larly polarised field with pulse envelope given by a general 
function [19|. The trident process was studied using the Weiz- 
sacker Williams approximation [20 1 and a complete calculation 
using pulsed laser field was performed by ll2Tl . One photon pair 
annihilation in a pulsed laser field with specific attention to op- 
erating parameters at the future XFEL was also studied recently 



The HICS process is also fundamentally important in ac- 
celerator physics where it is known as the Beamstrahlung and 
where it constitutes the main source of background photons due 
to beam-beam effects at the interaction point of a collider [23 1. 
The external field in the Beamstrahlung is provided by an on- 
coming relativistic charge bunch and is a constant crossed elec- 
tromagnetic field described by a 4-potential of infinite period. 
Interactions within a constant crossed field permit the radiation 
of vanishingly soft photons and the differential transition rate 
is divergent ll24l . 

In order to study the general properties of the HICS pro- 
cess it would be useful to have an analytic expression valid for 
all external fields for which exact solutions of the Dirac equa- 
tion are possible and of which those external fields described 
hitherto in this section are specific examples. The advantage of 
such a general analytic expression is that not only known ex- 
pressions for particular external fields could be written down 
immediately from a single source, but also the HICS transition 
probability for new external fields not yet studied theoretically 
would be readily available. The purpose of this paper, then, is 
to derive an expression for the HICS transition rate in a general 
external plane wave electromagnetic field and to illustrate its 
usage with particular examples. 



2 The Volkov Solution for a general plane 
wave external field 

In terms of notation, we employ a metric with signature 
(1,-1,-1,-1) and write the time and 3-vector part of a 
4-vector by x M = (a; , x). Conventionally, we work in natural 
units h = c = 1. 
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Fig. 1: The HICS process Feynman diagram. 



A general plane wave electromagnetic field with wavevec- 
tor k directed along the x 3 axis of a Cartesian coordinate sys- 
tem (figure |2l is described by a 4-potential A c ^ which can be 
decomposedinto two 4-vectors ain = (0, ai), a 2fl = (0,a 2 ) 
whose 3 vector parts are parallel to the xi,X2 axes and whose 
co-factors A?' ,A\' are general functions of arbitrary period 
described by a sum of N terms, 

N 

Al(k ■ x) = ai^Afik ■ x) + a^Af(k ■ x) (1) 



The strength of the external field is an important factor de- 
termining the extent of multi-photon interactions, and is char- 
acterised by the dimensionless parameter v which is a ratio (in 
natural units) of the fermion charge e, mass m and field poten- 
tial, 



en 
m 



a = ai = a 2 



(2) 



In the Furry picture the Dirac equation containing the 4- 
potential A^ for a fermion of mass m and momentum p = 

(e P ,p), 







(3) 



can be solved exactly, resulting in wavefunctions \P p first 
obtained by Volkov [25|. These Volkov solutions are a product 
of the normal free fermion solution including the bispinor u v 
and normalisation factor, with an extra phase S p and a term 
containing 4- e $, 



&p( k ■ x) = t= E p (k ■ x) u p 



(4) 



where E p (k ■ x) = ( 1 



exp I — i [ p + 

(k-x) r 



2(k ■ p) 

e 2 a 2 ( 



2(k-p) 
e(A e ((j})-p) e 2 A e (4>) 2 



x %Sp{k ' x^j 

e 2 a 2 ii 



(k-p) 



2(k-p) 2(k-p) 



In the Volkov solution above we have anticipated a 
quasi-momentum term containing the parameter £ which 
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has a linear dependence on the scalar product (k ■ x). The 
quasi-momentum is physically due to the average change in 
motion of the electron embedded in an external, oscillatory 
external electromagnetic field. The precise form of £ depends 
on the specific external field the HICS process takes place in. 




Fig. 2: The HICS coordinate system. 



The Volkov wavefunctions are used within the normal 
S -matrix theory to study interactions between the fermion 
embedded in the external field and the quantised boson field. 

In order to calculate a transition probability, the depen- 
dence on the scalar product (fc • x) must be simplified. The 
natural way of doing this is to write down a matrix element 
from a Feynman diagram (figure [T] with double straight lines 
representing Volkov wavefunctions) and gathering the depen- 
dence on space-time around the diagram vertices. The result 
is a modified vertex 7^ which can be transformed from posi- 
tion space to momentum space using a Fourier transform for a 
function of arbitrary period 2ttL, 

7£(P/>P*?fc- x ) = Ep f ( k ■ x hnE Pi (k ■ x) 



^(Pf,Pi,k-x)= / KZT ( 5 ) 

r=-oo •'— 

(ij; [4>v ~ (k ■ x)]^^(p f ,pi,4> v ) 



. exp 



3 The HICS transition probability for a 
general external field 

The HICS transition probability is calculated in the usual fash- 
ion, writing down a matrix element, squaring it, performing a 
trace and integrating over the phase space. 



fermion 4-momenta (£i,Pi), (e/,p/) and the radiated photon 
4-momentum (ojf,kf), the matrix element is 



M 



f w f 



E 



L 

-ttL 



d<j>v 
2tvL 



6 [pf + k f - pi 



e 2 a 2 £(k ■ kf) 
2(k- Pi )(k-p f ) 



k (6) 



Squaring the matrix element results in a double summation 
(over indicies r, r'). The subsequent product of delta functions 
will yield zero for r ^ r' and otherwise the product of space- 
time volume VT in which the process takes place and a single 
delta function in index r. The spin sums yield the usual pro- 
jection operators A p = f + m and trace, and there is a double 
integration over variables <p v , (j> w 

\M fl \ 2 _ 



E 



00 

E 



VT 32ne l e f u f ^ J (2nL) 2 

T — — OO 

iPf 1 v'idJJ 

. exp (i^r(4> v - (f> w )\ 
.( , ( r e 2 a 2 i(k ■ kf) \ , 

s i Pf + kf _ Pi _i - a(Jk . w)(Jb . w) )/, 



(7) 



The trace calculation is lengthy but straightforward. A sim- 
plification was obtained by grouping terms in order to separate 
the derivative with respect to the integration variables tfi v ,(j) w . 
The result of the trace is, 



4m 



(fc-Pi) , ( k -Pf) 
2m? \(k-p f ) (k-pi) 



[A*(<p v )-A e (<j> w )Y 



. exp it 



[S, 



Spi (4>v)] - i [S Pf (<f> w ) - S Pi (4> w )] ) 



(8) 



The integration over phase space is carried out as usual with 
the aid of the delta function. It is convenient to transform to in- 
tegrations over 4-vectors and perform the integration over the 
final electron momentum leaving the integration over the pho- 
ton 4-momentum, 



dkf dp f 



4e f lo 



$ [Pf + k f - pi 



e 2 a 2 j(k-k f ) \ 
2(k- Pi )(k-pf) 



= / d^kf 6{kj) 



.6 12 



e 2 a 2 Z(k ■ k f ) 
2(k- Pi )(k-pf) 



{k-p f )-2{ Pl -kf)\ (9) 



The matrix element for the HICS process can be written 
down directly from the Feynman diagram in figure [T] The ma- 
trix element is essentially the Fourier transform of the Volkov 
vertex with a Dirac delta function expressing the conservation 
of momenta, including a contribution j-k from the external 
field and the quasi-momentum term. Denoting initial and final 



In performing the integration over final photon momentum 
it proved useful to transform to light-cone coordinates, 



d 4 k 



f 



where k 



f 



-dki 
2 f 

-k^ 
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dk) dk) 



k° 



k 3 



(10) 
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Choosing the reference frame in which pi = for sim- 
plicity, the integration over dkj is performed using the sec- 
ond delta function in equation |9]). The integration over dkj is 
transformed to one over u = (k ■ kf)/(k ■ Pf) and is retained. 
The interim result for the transition probability is 



W 



In 



E 

T— — CO ' 
-i 2 1 ,2 „ 



w j, ^LJ k)dk2f 



7 



1,12 _ t.22 

K f K s 



e 2 l + (l + u) 2 



2m 2 



1 + u 



exp ( / 

</>, 



2(fc- Pi ) 



1/2,,, \ 



(fc r A e (0))^- 



(11) 



In simplifying the interim result above for the HICS transi- 
tion probability (equation ( fTTj )) there are two ways of proceed- 
ing. The first way, which proves useful for constant external 
fields, is to transform the delta function to an exponential and 
then to use the Fourier transform of a 2irL period function f(<p) 
to remove the dependence on r, 



E 



r— — oo 



ttL 



2ttL 



exp 



(i£[0-A])/W) = /(A) (12) 



The integrations over dkj dkj appear as Gaussian integrals 
and are readily carried out to obtain the HICS transition proba- 
bility which we label Wj^ s . Defining auxiliary functions, 



carry out the dkj dkj integrations in polar coordinates, zdz d9. 
The delta function from equation ( fTTj ) becomes, 



2 <z 



m 2 (l + u 2 £)u\ (k ■ kf)(k ■ Pf) 



2{k- Pi ) 



(15) 



and can be used to perform the integration over dz, thereby 
setting a range for the summation over r and restricting the up- 
per limit of the du integration. The integration over dd results 
in a Bessel function. It proves useful to make a change of inte- 
gration variables from dcp v d<j> w to d<j>- d(f>+. The end result is 

-2nd 
HICS 



labelled W&L and is 



2nd 
HICS 



E 



L 
-tsL 



du 



(2^L) 2 J Q (1 + u) 2 



2 + 1 + a + n)2 F 1 (cf> v , 0^) *> Jo ( 2zF 2 1 



2(1 + u) 



exp 



where, 



2(k- Pi ) 



e 20_ + « 



2(k- Pi ) 



2{k- Pi ) 



m 2 (l + u 2 £)u 
2{k- Pi ) 

(k.pi) r 



F, 



-1 1/2 



(16) 



+ 4>w) 



0(r) = for r < 0, <9(r) = 1 for r > 

We now have analytic representations for the HICS transi- 
tion probability for a plane-wave electromagnetic field of gen- 
eral form. It should be now possible to derive existing results 
from the literature and to write down new ones for external field 
forms not yet considered. 



Fi{ 



■0 



F : 



-[A%<j> v )-A^ W )Y 
A e ((l>)d(j) 
A e {(j)) 2 d(j) 



we have, 



w lst 

"HICS 



1- 



du 



T 2-kL L (1 + u) , 
ttL JQ v ' J— oo 



indX 
2X 



l + + 
2(1 + u) 



.exp i 



2(k- Pi ) 



F 2 (\ + 4>w,4>iv) 2 



(13) 



+ F 3 (\ + <j) w , <j) w ) 



(14) 



The second route from the HICS interim transition proba- 
bility, which proves useful for oscillatory external fields, is to 



4 The HICS transition probability for 
specific external fields 

To illustrate the general expression for the HICS transition 
probability obtained in the last section, two analytic results and 
a numerical result are obtained for particular external fields of 
practical interest. The specific analytic form for the HICS tran- 
sition probability can be obtained from either equation ( fT4| ) or 
( p"6*| ) after substitution of the particular 4-potential A e of pe- 
riod 2hujL into the auxiliary functions Fi,F 2 , F3, choosing a 
value for the quasi-momentum parameter £ in order to simplify 
the analytic expression, and by performing further integrations 
where possible. 



4.1 A circularly polarised external field 

Experiments investigating strong field effects are usually per- 
formed with intense lasers which commonly provide a circu- 
larly polarised electromagnetic field. Such a field is described 
by the 4-potential and resultant auxiliary functions, 
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A*((j>) = ai M cos(0) +a 2Al sin(0), L = 1, £ = 1 
Fx(<t> v ,<p w ) = 2v 2 (l -cos2</>_) 

F 2 ((f> v ,(f> w ) = 2z/|sin<£_| (17) 
Fz{4> v ,<j) w ) 



2v 2 d>- 



The quasi-momentum parameter £ = 1, results in a simpli- 
fied exponential dependence in the second general form of the 
HICS transition rate, equation (16) . There is no dependence in 
the transition probability integrand on the variable (j> + render- 
ing its integration trivial. The d(f>- integration can be carried 
out with the aid of the Bessel function identity (|26| §2.6), 



J r (z) 2 = 



2- 



Jo(2zsm0_) exp(i2r</>_) (18) 



The final result for the HICS transition probability for a 
circularly polarised external field is the same as that found in 
the literature (E3, equation 101.15), 



Tj/circ em x > 

Wmcs - -^r 2^ 



r=l 







2J 2 (z) - 



u 2 l + (l + u) 2 



du 
(l + «) 

2 



z = 2^ 



2(fc.pi) 



1 + M 

m 2 (l + i/ 2 )u 



(j2 +1 (z)+j2_ 1 (z)-2j2(z)) 

-,1/2 



2(fc-Pi) 



_ 2(fc. Pi )r 
~~ m 2 (l + u 2 ) 



(19) 



The above expression for the HICS transition probability 
has been extensively numerically studied in the literature. How- 
ever it is worth noting that a numerical study could proceed 
directly from the general expression for the HICS transition 
probability in equation ( |16| l. 



4.2 A constant crossed external field 

The Beamstrahlung results from charges within a particle beam 
at the interaction point of a collider decelerating and radiating 
within the strong field of an oncoming, ultrarelativistic charge 
bunch. Such an external field can be described by an infinite pe- 
riod, constant crossed electromagnetic field. The first form of 
the general HICS transition probability (equation ( fT4) ) will al- 
low us to directly write the known result for the HICS transition 
probability in this field. The 4-potential and resultant auxiliary 
functions are, 



Affl) = a^tj), L = oo, C = 
Fi{\+<f> w ,<f> w ) = v 2 X 2 

F 2 {\ + (t)w,4'w 



-A 2 



(20) 



F 3 (X+ip w ,<p u 



X A 



There is no dependence on the quasi-momentum parameter 
£, reflecting the fact that the electron does not oscillate in a 
constant field. Substitution of the infinite period L = oo, poses 
no problem since there is no dependence of the integrand on 
the d(j) w integration and, 



-7r£ 



(2ttL) 



(21) 



Substitution of the auxiliary functions of equation pO] ) 
yields an interim HICS transition probability for a constant 
crossed external field, 



Theorist _ 
"HICS — 



■Itt 



i dX 



du 

o (i + «) 2 

1 v 2 l + (l + u) 2 



A 



, exp i 



1 + u 



2(k 



■Pi) V 12 J J 



(22) 



We need now only to rescale A to get the integral repre- 
sentation with respect to A of an Airy function. The final form 
of the HICS tranition probability for a constant crossed field, 
which is the same as that appearing in the literature ESI , is 



TTAConst 
K|/ H1CS 



(i + «) 2 



dz ■ 



1 + (1 + m) 2 d 



z{l + u) dz 



where 



v{k ■ Pi) 



2/3 



Ai(z) 
(23) 



4.3 Two phase offset circularly polarised external 
fields 



In at least two practical endeavours, namely laser-plasma ac- 
celeration of charged particle bunches [29 1 and laser confine- 
ment fusion experiments l30l . multiple intense laser interac- 
tions with matter take place. In that spirit we examine the HICS 
transition probability for the simple multi-field case of two co- 
linear circularly polarised external fields of equal energy and 
intensity, separated by a variable phase difference denoted by 
X- The 4-potential of the combination of field components and 
the resultant auxiliary functions are, 

Al{4>) = ai M (cos0 + cos(0+x)) + a 2u (sm^) + sin(0+x)) 



L = l, e = 2(l + cos X ) 

4>w) 



2 

,4> w ) = Av cos 

2 



v 2 cos 2 { ? ) (1 - cos2</>_) 



(24) 



sin< 



F 3 ((t>v,<t>w) = 4^ 2 (1 + cosx)4>- 

The quasi-momentum parameter £ in this case depends on 
the phase difference between the two fields. After substitution 
of the auxiliary functions into the second general expression 
for the HICS transition probability a transition probability for 
two phase-offset circularly polarised fields, W2° cs is obtained. 
Apart from the appearance of the phase difference \ me an_ 
alytic form of the transition probability resembles that for a 
single circularly polarised field, 
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where 



z = Av cos : 



■ (j; + ,w+j;-,w-2j;w) 



m 2 (l + 2;/ 2 (l + cosx))m 
2(fe- Pi ) 



(25) 



1/2 



2(fc- Pi ) 

2(fc.pQr 
m 2 (l + 2^ 2 (1 + cosx)) 



A numerical evaluation of equation ( |25j l is straightforward 
since the summand converges quickly and the oscillatory na- 
ture of the integrand is inhibited. Figure [3] shows the variation 
of the HICS transition probability with phase offset for differ- 
ent external field intensities. The HICS transition probability is 
maximal when the two fields are in phase, and vanishes when 
the phase difference is tt and the two fields destructively inter- 
fere 




co-linear components have been used. By choosing such a 
representation for a general plane wave, any waveform that 
can be expressed by a Fourier series can be studied using 
general forms of the HICS transition probability written down 
in equations ( fT4| i and ( p"6| ). 

The validity of the general HICS transition probability 
expressions was established by substituting in particular 
4-potentials - that of a circularly polarised and of a constant 
crossed electromagnetic field - to obtain the same analytic 
forms of the HICS transition probability found in the literature. 
Specific numerical studies can, however, be launched directly 
from the general HICS transition probability expressions 
which is particularly useful for external fields for which no 
simple analytic expression for the transition probability is 
possible. 

A multi-external field case, in which the HICS process 
takes place within two co-linear circularly polarised fields of 
equal intensity and energy separated by a phase difference, 
was examined analytically and numerically. As the phase 
difference reaches tt, the two fields destructively interfere and 
the HICS transition probability vanishes. 

New high intensity laser facilities are being planned and 
built which will provide the means to experimentally test 
physics processes in strong electromagnetic fields. A variety 
of experimental conditions involving different configurations 
and combinations of external fields will be possible. In the case 
of experimental tests of the HICS process we have provided 
a general expression for the HICS transition probability from 
which, no matter what the experimental configuration, specific 
analytic expressions and numerical values can be quickly gen- 
erated. 
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